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We investigate a lattice of interacting anharmonic oscillators by using a mean field theory 
and exact diagonalization. We construct an effective five-state hopping model with intersite 
repulsions as a model for /3-pyrochlore AOs206(A=K, Rb or Cs). We obtain the first order 
phase transition line from large to small oscillation amplitude phases as temperature decreases. 
We also discuss the possibility of a phase with local electric polarizations. Our theory can 
explain the origin of the mysterious first order transition in KOS2O6. 
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Recently, various unusual low-temperature properties 
of low-energy anharmonic phonons have attracted much 
attention. In /3-pyrochlore compounds AOS2O6 (A=K, 
Rb or Cs),^"^ low-energy anharmonic oscillations of A- 
ions are considered to play an important role for the real- 
ization of the unusual electrical resistivity-'^'^ and specific 
heat,^ the strong coupling superconductivity,^"^ and the 
nuclear magnetic relaxation time,^ especially in KOS2O6. 
In these compounds, A-ions are located in large OS12O18 
cages and form a diamond lattice structure. The cage 
opens in the four [111] directions, which are the bond 
directions of the diamond lattice, and it is expected that 
A-ion oscillations are anharmonic and have large am- 
plitudes along these four directions.^' ^ Since the atomic 
radius of K is the smallest among A-ions, the effects of 
the anharmonicity is most prominent in KOS2O6. 

One of the most interesting phenomena in KOS2O6 
is the first order transition below the superconducting 
transition temperature Tc — 9.6 K.^ The transition tem- 
perature Tp =7.5 K is robust upon applying magnetic 
fields up to 14 T.^^ This suggests that the transition is 
driven by lattice degrees of freedom. Interestingly, it is 
found that, in the non-superconducting state under mag- 
netic fields, the temperature (T) dependence of resistiv- 
ity changes from an unusual \/r behavior to a depen- 
dence expected in the fermi liquid theory.^ Correspond- 
ing to this, the life time of quasiparticles becomes longer 
below Tp}^'^^ The energy gap of the superconductivity 
observed in the photoemission spectroscopy also shows 
an anomaly at Tp}"^ These experimental results suggest 
that quasiparticles are coupled to the lattice degrees of 
freedom leading to the first order transition. 

As for the theories, it is proposed from the band struc- 
ture calculation, that the effective potential of A-ion is 
shallow and even has a minimum at off-center positions, 
especially in KOs206.^^ Dahm and Ueda investigated the 
origin of the unusual temperature dependence of the re- 
sistivity and the saturation behavior of the NMR relax- 
ation rate at high temperatures.^'* As for the first order 
transition, Kunes et al, proposed that it is multiple-q or- 
dering of the potassium displacements.*^ This scenario, 
however, contradicts the observations that symmetry is 
not broken below Tp as various experiments show."^' ^' 

In this Letter, we extend the model used in Ref. 15 to 



include quantum hoppings by introducing five localized 
states of the ionic configurations. One of the five states is 
an on-center state and others are off-center ones. We an- 
alyze this model by using a mean field theory and exact 
diagonalization (ED) calculations. In our theory, the first 
order transition is characterized by the A-ion density at 
the on-center position, which is an order parameter in 
a sense of a liquid-gas transition. We will discuss impli- 
cations of the present results for KOS2O6 and physical 
consequences of this phase transition. 

Before going into detailed calculations, we discuss 
the first order transition without symmetry breaking in 
K0s20g from a phenomenological point of view. The lo- 
cal symmetry around A-ion is point group Td and lo- 
cal potential V{r) of A-ion is approximated as V{r) — 
a|r|2 -\- bxyz -h c|r|^ -\- c'{x^ + + z") + 0(|r|^), where 
r = (x,y,z) is the coordinates of A-ion from the on- 
center position of the diamond lattice. It is noted that the 
third order term bxyz can exist. Due to this term, V{r) 
becomes shallow in the four [111] directions of the bonds 
in the diamond lattice. This makes s-wave (ips) and f- 
wave (ipxyz) functions of oscillation states hybridized. 
Both of them belong to Ai representation in point 
group. The one-body ground state ^E" for the potential 
V(r) is a linear combination of tps, fxyz and other Ai 
components: vE" = Csfs + Cxyz^xyz + • • ■ • The point is 
that ^' has inner degrees of freedom characterized by 
the coefficients Cs and Cxyz, etc. These coefficients can 
vary without breaking any symmetry. In this Letter, we 
propose a minimal model to exhibit such variation as 
functions of T and flatness of V{r). When this variation 
becomes a sudden jump, the first order transition occurs. 
In this scenario, the 6-term plays a more important role 
for anharmonicity and anisotropy than the fourth order 
terms. It is noted that this kind of transition is com- 
pletely different in nature from classical structural tran- 
sitions, which are described by the instability of some 
phonon modes. 

Now, we introduce a minimal model to describe the 
phase transition at Tp in KOs206- We investigate a 
Hamiltonian of five localized states with inter-site re- 
pulsive interactions which are originated from Coulomb 
interactions between A+ ions. The five localized states 
are one on-center = 0) and four off-center = 1,2,3 
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Fig. 1. (a) Cubic unit cell of diamond lattice with one on-center 
and four off-center sites. On- and off-center positions are indi- 
cated by large and small filled circles, respectively, (b) Interac- 
tions and tunneling matrix elements. 

and 4) states as shown in Fig. 1. These five states can 
describe two kinds of Ai components: (ps and ipxyz and 
also p-wave hke components (ppi{i — x,y and z) which 
have the same symmetry T2 as one-phonon states. The 
on-center state is located exactly at the diamond lat- 
tice site — {xia,yia, Zia) wherc i belonging to sublat- 
tice a{= A or A = B). Others are located at positions 
XjQ + shifted along the bond directions of the dia- 
mond lattice. We can regard that the parameter \da(\ is 
implicitly included in parameters in our model Hamilto- 
nian shown below. We note that 6-term in V{y) prefers 
wavefunctions having large weight in the bond directions, 
and thus Coulomb interactions for these wavefunctions 
become effectively large. These bases rather than phonon 
ones are expected to be a good starting point to discuss 
AOS2O6, since the local oscillations of A- ion is strongly 
anharmonic and large in the directions of the bonds es- 
pecially for A=K. This choice of the bases corresponds 
to write 4" in the form of Ca'^a + Ca/'^^a' ■ Here, (fa and Lpa' 
correspond to the on-center state and a symmetric linear 
combination of off-center ones, respectively. 
The Hamiltonian in this Letter is written as 



H 



(1) 



where 0?^^ represents a creation operator of the ionic 
state jia^) localized at position ^ at a site i of sublattice 
a and Uia^ = ^J^^^ia^. The five states labeled by ^ and 
rj have hopping matrix elements e^^ with ^ rj which 
represent the local quantum kinetic energy. We set = 
-t for ^ > 1 and e^^ = -t' for ^,7? > 1. For the local 
energy levels, we set Eqq — and e^^ = e for ^ > 1. For 
the interaction terms, we use only the nearest neighbor 
Coulomb repulsion U and the next nearest one < C/' < 
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Table I. Eigenstates and eigenenergies of the mean field Hamil 



tonian. e 
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'5)/v'(^ - ^ - ^00)^ + 4*^- The wavefunctions are written as 
ipa = <t>0, Va' = \{4>l+<l>2+<l>3+<t>'i), ^Px = 5 ((/'I - </'2 + 03 " </'4 ) , 
ipy = \{li>\ 

where we omit the site indices. 



/>2 - (/is + 04), and = 5 (01 -I- 02 - 03 - 04), 



- U/t=6, UVt=5 






t'/t= 


D.2 1st...<;>* 


2nd 




\ y'^ st 


"2nd «>,^ 




SBP / 


/I St 


^ t'/t=0.4 









-1 







8/t 



1 



Fig. 2. Mean field e-T phase diagram of model (I). U = 6t, U' = 
5t and t' /t = 0.4 and 0.2. 

U as shown in Fig. 1 (b). This is because the screening 
by the electrons on the cage is not effective for these two 
terms because of the peculiar structure of the cage, but 
effective for others, and thus we neglect them. 

First, we investigate this model by a mean field approx- 
imation (MFA) in the case that the translational symme- 
try is not broken. The mean field decoupling of the inter- 
actions leads the effective local energies Cqq = AU'{na) 
and 6°" — e + U (ng) -|- U'{nao), if the local symmetry 
is not broken. Here, (•) denotes the thermal average and 
(ria) — (riia^) with ^ > 1 and (n^o) — ("-iao)- Now we 
can easily diagonalize the mean-field Hamiltonian, and 
show the eigenenergies and the eigenstates in Table I. 
There are two singlet (Ai) and one triplet (T2) states. 
The ground state A^^"* is an in-phase combination of (pa 
and (fa' but the excited A\ state is an out-of-phase one. 
Using these eigenstates, we calculate (nao) = 1 — 4(nQ) 
such that the self-consistency is satisfied. 

When we consider a symmetric phase with (n^g) = 
(riBo) = no at T = 0, the self-consistency condition reads 



(2) 



with/3(no) - V(l - ^o)/no and /(no) = [e + U/4-U'- 
3t' + {2U' ~ C//4)no]/4t. A first order transition at T = 
occurs if U' — U/8 > At as e is varied, since this condition 
is satisfied, Eq. (2) has three solutions. By using Eq. (2), 
the ground state energy per site Eg is represented by the 
function of no as Eg = (U' - C//8)(l - no)^ - 2tp{no). 

In Fig. 2, we show the mean field e-T phase diagram 
for U/t = 0.6, U' = 0.5t and t' /t = 0.4 and 0.2. There are 
the first order transition line, and the second order tran- 
sition to a symmetry broken phase (SBP). Between the 
two transitions, there is a region where no remains ^1/5 
down to r = for t'/i = 0.4. In this region, the oscilla- 
tion of A-ion becomes large and gains the kinetic energy 
rather than the potential energy. This intermediate re- 
gion becomes smaller as t' decreases and eventually the 
FE phase touches the first order transition line. When 
we increases U, the former is suppressed as understood 
from the T = analysis above and becomes a crossover. 
On the other hand, the transition temperature to SBP 
is enhanced. This is because only U enters in the self- 
consistent equation for SBP, and the other parameters 
are implicitly included in it through the energy levels. 
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Fig. 3. Temperature dependence of (a) ng and (b) entropy. U = 
6t, U' = 5t and t' /t = 0.4. 



First, we discuss the first order transition which is rel- 
evant to /3-pyrochlore compounds, and then investigate 
the properties of SBP. Note that we consider the case of 
e < is not reahzed in AOS2O6. 

In Fig. 3 (a), the temperature dependence of no is 
shown for various values of e. At high temperatures, the 
density no is about 1/5. This means that the five posi- 
tions are equally populated due to thermal fluctuations. 

As for the energy level, for e = 0.8t, ct^ ^ ^ aW ^ 3i 

^. 

and e .(2) — e «(i) ^ 5< in the high temperature limit. In 
the case of large e, there is no discontinuous change but a 
smooth crossover to a large no value with decreasing tem- 
perature. As e decreases, the crossover becomes a jump, 
i.e., the first order transition. Just below the first or- 
der transition temperature T*, no suddenly increases in 
order to reduce the potential energy e and the intersite 
repulsions. Consequently, ejvj increases below T* . This 
means that the effective potential around the central po- 
sition becomes steeper below T* . This is understood by 
considering the T ~ Q case. When U' — U/8 > 0, which 
is automatically satisfied by that for the first order tran- 
sition at r = 0, increases as no increases. 

In Fig. 3 (b), we show the temperature dependence of 
the entropy. The entropy at high temperature approaches 
to log 5 ~ 1.6. There is a jump of entropy for 0.6 < e/t < 
1.4, accompanying the first order transition. The released 
entropy at T* is ^ 5 log 2 per A-ion for e/t = 1.0. The 
jump of entropy becomes small as e increases from e/< = 
1.0 and disappears at the critical end point e ^ I. At and 
T ~ 1.6i and also at T* 0. 

In the SBP, the average {niAi + n\A2 — "-i/is — niAA) = 
{niBi+niBi — n-iBi — niB2) becomes finite uniformly. This 
means a finite polarization P || [001], i.e. ferro-electric 
(FE) state. This phase is stabilized when T2 states have 
a low excitation energy. This corresponds to small and/or 
negative e. The reentrant behavior as a function of tem- 
perature is due to the fact that the FE fiuctuations arise 
from the excited T2 states. It is noted that the FE state is 
six-fold degenerate (P || ±[100], ±[010] and ±[001]). For 
the three-dimensional six-state Potts model, it is known 
that the transition is the first order both from MFA and 
Monte Carlo simulations.^^ Our results of MFA is differ- 
ent from the result of MFA for the Potts model. This is 
because the symmetry P ^ — P is present in this model. 

We have examined all the possible uniform symmetry 
broken patterns in our MFA. The phase with uniform po- 
larizations is the most stable one. For example, a phase 



where the i?-sublattice polarization is antiparallel to the 
A-sublattice one is not realized, because of repulsive in- 
teractions. The energy of another FE phase with P || 
[111] is slightly higher than that of P || [001]. When we 
consider an order with the wavevector Qj, ~ (27r, 0,0), 
the order parameter is {niAi + niA2 — 'T-iAs ~ n-iAi) — 
{niB2 + n-iBs — riiBi — ^iS4) 7^ 0. This corresponds to 
polarization P^ || [001] and P_b [| [010] in the two sub- 
lattices. This state is one of the twelve-fold degenerate 
states; P^ —Pa, A B and Q^; — > Qj, or Q^. The 
energy of these states become lower than that of the FE 
state when long-range repulsions are included. 

Next, in order to investigate the effect of ffuctuations 
beyond the MFA, we carry out ED calculations of this 
model. We use the cubic unit cell which includes eight A- 
ions and retains the Td symmetry with periodic boundary 
conditions. Since there is no conserving quantities such 
as spin or charge as in electron systems, the ED in larger 
system size is difficult at present. However, we can see 
the effects of the fiuctuations in this size. 

The expectation value ni^o = ^Og in the ground state 
is shown in Fig. 4 (a) as functions of U and e. Here, nog is 
independent on the site indices. We find a sharp change 
of nog for relatively larger U than in the case of mean field 
results. This may be due to the effects of the quantum 
fiuctuations and finite size. In this calculation, the sharp 
change of nog (e) is due to the anti-crossing of the ground 
and the first excited state as shown in Fig. 4 (b). Since 
the off-diagonal elements of the Hamiltonian (1) are all 
non-positive, we can apply the Perron-Frobenius theo- 
rem and prove that the ground state is always unique. 
This means there is no level crossing for the ground state. 
We expect this anti-crossing becomes narrower with in- 
creasing the system size, and evolves into a first order 
transition in the infinite size limit. 

The presence of SBP is also reficcted in the structure 
of excited states. When e is smaller than ~ t, the first 
excited states are nine-fold degenerate as shown in Fig. 
4 (b). As discussed before, six Q = states and twelve 
Q ^ states are degenerate in the MFA. In this finite 
size cluster, these eighteen states are coupled and split 
into a few multiplets and the lowest excited multiplet has 
nine states. This unusual large degeneracy is due to the 
special topology of eight-ion cluster. Each A-sublattice 
site is connected to all _B-sublattice sites, and vise versa, 
and the total symmetry turns out to be ^ Z2, 
where and Z2 represent iS'4 group for sublattice a 
and the permutation A ^ B, respectively. This special 




Fig. 4. (a) Occupation number nog as functions of U and e. (b) 
Energy eigenvalues of four low-energy multiplets. Numbers rep- 
resent the degeneracy. U/U' = 1.2 and t' /t = 0.4. 
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symmetry is broken in larger system sizes and it remains 
as an interesting question how the degeneracy is hfted. 

The tendency of the first order transition from small to 
large riog a-s e increases, is consistent with the results of 
MFA and due to the quantum fluctuations, the transition 
occurs at larger U and U' than in the MFA. In order to 
conclude whether the intermediate region obtained in the 
MFA, i.e., the region e for t' /t = 0.4 in Fig. 2 exists 
or not, it is needed to carry out detailed analysis of the 
finite size scaling for the energy of excited states. 

As we discussed above, this model shows the first or- 
der transition at T* below which the on-center density 
suddenly increases. We believe that this corresponds to 
the first order transition at Tp in KOsaOe- Since anhar- 
monicity for Rb and Cs compounds is weaker than that 
for K, we can regard that the cases for Rb and Cs corre- 
spond to the larger e region in the phase diagram where 
there is no transition but a smooth crossover. The fact 
that this transition does not break any symmetries is 
consistent with experimental results.^' ^'^^ Recently, Ya- 
maura et al., observed that the charge distributions of 
potassium shrink below Tp by the neutron experiments 
in KOsaOe-^^ This corresponds to the increase of uq. 
In addition to this, the suppression of the quasiparticle 
scattering below Tp^' is explained as a result of the 
sudden freezing of the A-ion oscillations. It is also noted 
that T* is typically one third of for e = 0.8t in Fig. 2. 
This is qualitatively in good agreement with estimated 
Einstein energy ~ 22 and Tp = 7.5 K. 

In real materials, the cage is composed of O and Os 
ions. Dynamics of these ions will be affected by the tran- 
sition at Tp. In /3-pyrochlore structure, ions can 
move without breaking lattice symmetry. Their positions 
are expected to be different in high- and low-temperature 
phases. When we take into account the effect of the near- 
est neighbor K+-ion from an 0^~-ion, the is ex- 
pected to move toward the K+-ion below Tp due to the 
increase of the attractive Coulomb energy. 

An interesting prediction for KOS2O6 of our theory is 
that the Tp increases when we apply pressme. This is 
understood by noting that the energy level of off-center 
states e is expected to increase, since e corresponds to 
the steepness of the potential V{r). It is interesting to 
explore the critical end point of the first order transition 
at high pressure. The critical exponents of this transition 
are expected to be those of the three dimensional Ising 
model. Interestingly, It is observed that the pressure de- 
pendence of the superconducting transition temperature 
Tc has a small anomaly at aroimd 2 GPa and T=8 K.-^^ 
This might be a signature of enhanced Tp. 

A recent theoretical study of the strong coupling super- 
conductivity for KOS2O6 pre- assumes the energy of the 
low-energy Einstein phonon decreases below Tp.^° This 
assumption is not consistent with our results, since the 
energy level of the T2-states increases below T*. The tem- 
perature dependence of the energy gap of the supercon- 
ductivity across Tp^^ might be explained by the balance 
between the increase of the phonon energy and the re- 
duction of coupling constants with conduction electrons. 

Theoretically, it is important to extend the present 



model to that can describe the higher energy physics, 
including the couplings with conduction electrons^^ and 
acoustic phonons. The validity of using the short range 
interaction should also be explained by a microscopic 
theory, since if one uses the bare Coulomb or a simple 
Yukawa type interaction, the first order transition does 
not occur in the MFA discussed in this Letter. Under- 
standing the first order transition in the phonon basis 
and the observed unusual T^ specific heat in KOs206'^ 
still remains as a big challenge and needs further theo- 
retical studies. 

In conclusion, we have investigated the origin of the 
first order transition in KOS2O6. We have proposed that 
the transition at Tp in KOS2O6 is an on- and off-center 
transition without symmetry breaking which is liquid-gas 
type. Our results show that such a transition is possible 
in the case of rather shallow potentials. Since the poten- 
tial in KOS2O6 is expected to be the shallowest in three 
compounds of AOS2O6 (A=K, Rb and Cs), this explains 
why only KOS2O6 shows the first order transition. 
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